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PARTICLE FLOW WITHIN A TRANSONIC COMPRESSOR ROTOR 


PASSAGE WITH APPLICATION TO LASER-DOPPLER VELOCIMETRY 

by Barry R. Maxwell 
Bucknell University 


SUMMARY 

A theoretical analysis was conducted of the dynamic behavior of 
micron size particles moving in the three-dimensional flow field 
of a rotating transonic axial-flow air compressor rotor. The 
particle velocity lag and angular deviation relative to the gas were 
determined as functions of particle diameter, mass density, and 
radial position. Particle size and density were varied over ranges 
selected to correspond to typical laser-Doppler velocimeter (LDV) 
flow field mapping applications. It was found that the particles 
move essentially on gas stream surfaces and that particle tracking 
is relatively insensitive to the rotor radial coordinate. Velocity 
lag and angular deviation increased whenever particle size or mass 
density increased, and particle tracking was more sensitive to a 
change in particle diameter than to a corresponding change in mass 
density. Results indicated that velocity and angular deviations 
generally less than 1 percent and 1 degree could be achieved with 
1 gm/cc tracer particles with diameters of 1 micron or less. 





INTRODUCTION 


Laser-Doppler Velocimetry (LDV) is a noncontact optical 
technique for determining the velocity of flowing fluids by meas- 
uring the velocity of tracer particles entrained in the fluid 
rather than by measuring the velocity of the fluid itself. The 
LDV technique is of considerable importance in gas dynamic meas- 
urements because it is unnecessary to introduce any hardware into 
the flow system. Because of the absence of a physical probe the 
LDV technique is ideally suited to measuring gas velocities in the 
flow passages between rotor and stator blades of turbomachines. 
Development of noncontact velocimetry techniques in turbomachinery 
will provide data currently required to, extend and improve design 
methods , blade performance and efficiency of components used in 


advanced airbreathing engines. It is for this reason that the LDV 
technique is beginning to find application in the aircraft turbo- 
machinery field (ref. 1,2) . 


The laser-Doppler velocimeter system measures only the motion 
of either naturally occurring or artificially generated micron size 

tracer particles entrained in the fluid i Its accuracy is, therefore , 

( 

limited by the accuracy with which the pjarticles follow the fluid 

i . 

I 

flow. Since it is important to know hov< accurately the particles 


track the fluid, careful consideration mjist be given to the magnitude 

• • ' ' J 

' . • • ' , j 

and direction of the particles velocity (relative to the carrier fluid 
velocity. This report documents a theoretical investigation of the 
three dimensional dynamic behavior of micron size particles entrained 
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in gas flowing through a rotating transonic axial-flow air com- 
pressor rotor. The particle-to-gas velocity ratio and angular 
deviation relative to the gas are determined as functions of 
particle diameter, mass density, and radial position. Information 
of this nature is currently required for the application of LDV 
techniques to intrablade flow field mapping of aircraft turbo- 
machinery. 

Several previous studies (ref. 3,4,5) have investigated par- 
ticle tracking for LDV applications but have been restricted to 
particle flow in such systems as turbulent flows, nozzles, oscil- 
lating gases, shock waves, and over airfoils. Reference 6 used a 
LDV system to measure the velocity lag of aerosol particles in 
small supersonic nozzles and then used this information in conjunc- 
tion with numerical predictions to determine particle size. A 
theoretical study by Hussein and Tabakoff (ref. 7) and an exper- 
imental study by Tabakoff et al. (ref. 8) provide trajectories 
and velocities of solid particles suspended in a gas passing through 
an axial flow turbine stage. The objective of their investigations 
was to estimate the flow properties of a gas-particle suspension 
flowing over turbine blades. This information, would ultimately be 
used to minimize blade erosion in dusty environments through appro- 
priate redesign. Since their work considered particle diameters 
generally ranging from 100-900 microns, their results are inappli- 
cable for LDV applications because of the resulting high velocity 
lags. Wisler and Mossey Cref. 2) obtained detailed flow field 
measurements within the rotating blade row of a low speed axial 
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compressor using a laser-Doppler velocimeter. Their measurements 
generally agreed to within 2-3 percent with hot film and pitot- 
static probe measurements outside the blade row, and with potential 
flow solutions near the leading edge. Their work was restricted 
to flow velocities sufficiently low that particle velocity lag was 
not regarded as a significant problem. However, as gas flow veloc- 
ities and accelerations increase, particle velocity lag in blade 
channels becomes more significant, and application of LDV techniques 
to flow field mapping should be guided by the results of a dynamic 
analysis of similar gas-particle flow systems. 

A two-dimensional dynamic analysis of micron size particles 
entrained in high Mach number subsonic gas flow passing through 
a stationary cascade of turbine stator blades was conducted by the 
author in reference 9. The particle-to-gas velocity ratio and 
particle angular deviation relative to the gas were determined as 
functions of particle size and mass density. Results indicated 
that velocity and angular deviations generally less than 1 percent 
and 0.5 degree respectively, could be achieved with 1 gm/cc tracer 
particles with diameters of 0.5 microns or less. The results were 
restricted, however, to a nonrotating blade configuration with a 
constant meridional streamline radius, and were not necessarily 
applicable to rotating systems with variable streamline radius. 

GAS -PARTICLE FLOW SYSTEM 

The flow system analyzed in this investigation is the three- 
dimensional gas-particle flow through a transonic axial-flow air 
compressor rotor. The rotor is currently being tested at the 
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NASA-Lewis Research center as part of a program v to reduce the size 
and weight of fans and compressors while maintaining high levels 
of performance. The rotor was designed for a weight flow per unit 
annulus area of 198 kg/sec-m 2 (40.6 lbm/sec-f t 2 ) , a tip speed 
of 424 m/sec (1393 ft/sec) , a total pressure and temperature ratio 
of 1.601 and 1.162 respectively, and a tip solidity of 1.30. This 
resulted in a rotor with '44 blades each with a 2.5 aspect ratio, 
a leading edge tip diameter of .5038m (1.653 ft.), a hub-to-tip 
radius ratio of 0.504, and a rotational speed of 16,100 revolutions 
per minute. A photograph of the rotor is shown in figure 1. Per- 
formance data ,and blade-element design parameters are presented in 
reference 10. 

Analysis of the gas-particle flow through the transonic com- 
pressor rotor in figure 1 is conducted from a radius of 60 to 94 
percent of the blade span as measured from the rotor axis. Analysis 
was not made below 60 percent span because each blade has a vibration 
damper located at about 55 percent span which precludes the use of 
LDV techniques below that level. 

TECHNICAL APPROACH AND ANALYSIS 

Analysis of the gas-particle flow field was conducted in two 
separate steps. In the first, the magnitude and direction of the 
gas velocity was determined at all points of a blade channel in 
the 60-94 percent span range. In the second, the magnitude and 
direction of the particle velocity was determined throughout the 
region by analyzing the path followed by particles in a three- 
dimensional flow field with a known gas velocity distribution. 
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Decoupling the gas and particle flows in this way is justified 
because of the low mass fraction of naturally occurring or arti- 
ficially generated aerosols in most LDV gas-dynamic applications. 
Therefore, since the particles do not influence the gas flow, the 
particles and gas may be analyzed separately. 

GAS FLOW ANALYSIS 

The three-dimensional gas flow field was determined with a 
numerical technique developed by Kurzrock and Novick (ref. 11) 
which yields a solution for steady, transonic blade element flow 
in an axial compressor rotor. The equations governing the flow in 
a rotating blade row are formulated as an abbreviated form of the 
Navier-Stokes equations with the dominant viscous terms retained 
to allow for the formation of embedded shock waves in the flow field. 
The coordinate system used is an orthogonal curvilinear system 
attached to the rotor. The independent variables are the meridional 
streamline distance m and the tangential angle 9 . For purposes of 
analysis, finite solution regions are chosen on b lade- t'o-b lade 
surfaces. These regions are bounded on the top and bottom by adja- 
cent blade surfaces and by arbitrarily chosen upstream and down- 
stream boundaries. Figure 2 illustrates a blade-to-blade surface 
of revolution, a typical solution region, and the coordinate systems 
used in this study. A stream channel is defined by specifying a 
meridional streamline radius and the stream channel thickness as 
functions of the meridional distance. It is assumed that the 
radius varies linearly with the axial distance. Therefore, the 
angle between any stream surface and the rotor axis as measured 
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in the meridional plane is constant and- is referred to as the cone 
angle X . 

Gas flow field analysis was conducted on stream surfaces at 
60, 69.2 , 77.8, 86.2, and 94.3 percent radial span. The governing 
equations are derived in reference 11 for the meridional-tangential 
relative frame of reference shown in figure 2. A thermally and 
calorically perfect gas has been assumed with constant viscosity 
and thermal conductivity. Shock waves are handled as a viscous 
phenomenon and, therefore, become "smeared" in the solution. The 
method of characteristics is used to determine the boundary condi- 
tions along the blade surfaces and at the exit plane. Flow tangency 
is enforced on the blade surfaces, and a constant back pressure, is 
specified at the exit plane. Steady supersonic uniform flow condi- 
tions are specified along the upstream boundary, and periodicity 
is enforced on the remaining free surfaces of the solution region. 

The entire numerical technique is incorporated into a computer 
code designated KOTOR. The code is written in FORTRAN IV for use 
on the IBM 370/165 and CDC 6600 digital computers, and was modified 
in this study for use on the Xerox Corporation Sigma 7 computer. 

The input to ROTOR requires the m- 9 blade coordinates on the stream 
surface being analyzed, and the compressor flow conditions. The 
output includes the relative meridional velocity, tangential velocity ^ 
flow angle measured from the axial— radial plane, and Mach number 
at each grid point of a network established for finite difference 
analysis. 
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PARTICLE FLOW ANALYSIS 


From Newton's second law of mechanics the vector equation of 
motion of a single, spherical particle moving with an absolute 
velocity Ug is 

ID 3 Pp 3p = P g V (Ug-Up)? (1) 

The term on the left denotes the mass times the absolute acceleration 
of the particle and represents the force required to accelerate 
the particle in the absence of potential force fields. The term • 
on the right represents the total viscous drag exerted on the par- 
ticle by the surrounding fluid and is expressed in terms of an 
empirically determined viscous drag coefficient C D . Equation (1) 
assumes that the force on the particle caused by the pressure 
gradient in the surrounding fluid, the force required to accelerate 
the apparent mass of the particle relative to the fluid, and the 
force which accounts for the deviation of the flow from steady 
state (Basset force) are negligible. These assumptions are justi- 
fied (ref. 5,7) in flow systems whose properites are changing very 
gradually or when the particle density is much larger than the 
fluid density. Both of these conditions are generally satisfied 
in this study. It is also assumed that the particles are non- 
interacting and uniformly distributed within the carrier gas and 
that there are sufficiently few particles that their presence does 
not alter the gas properties from that of the gas flow alone. 

These assumptions are valid for most LDV applications employing 
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naturally occurring or artifically generated aersols. 

Letting C D - 2 4^(Re, M) /Re, where \|>(Re, M) is a Reynolds-Mach 
number dependent function which accounts for deviation of the drag 
coefficient from the ideal Stokes flow value (24/Re) , equation 
(1) becomes; 

a p = e <“g- " P ) (2) 

where 3 = 18vn]j(Re >M) . 

D 2 p p 

In order to derive the differential equations of motion of a 
particle passing through a rotating blade row of a compressor 
rotor, the coordinate systems o'lr 0 and oxyz are defined and are 
illustrated in figure 2, The frame of reference o'lr 0 is fixed 
in an arbitrary blade passage of the rotor and moves with constant 
angular velocity w of the rotor. The coordinate directions 1, r, 
and 9 are the axial/radial, and tangential directions respectively, 
and are used to denote the location of a particle relative to the 
rotating blade row. Point O' is the origin of the relative coor- 
dinate system and is located on the rotor axis at the intersection 
with a perpendicular plane passing through the leading edge of the 
solution region. The oxyz coordinate system is an inertial frame 
of reference , attached to the rotor axis and is composed of mutually 
perpendicular axes y, and z. The x coordinate axis is in the 
axial direction (1 direction of the relative coordinate system) 
while axes y and z lie in a plane parallel to the r - 0 plane of 
the relative system. 
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Assume that at time t=o the planes formed by the 1-r axes 
and x-y axes are coincident and a particle p is at an arbitrary 
position 1, r, and 0 as measured from o’. As illustrated in 
figure 2, after time t the o’lr system has rotated about the 1 
axis, and the 1-r plane is an angle of u>t from the x-y plane. 

At any point in time the x, y, and- z components of the particles 
equation of motion are, from equation 2: 


a px " & ( u gx “ u px^ 
a py = 3 (ugy - Upy) 

apz •= 3 (ugz - Up Z ) — ^ 

The relationships between the absolute and relative axial compo- 
nents of the particle displacement, velocity and acceleration are 
seen from figure 2 to be: 


Xp — 1 

u P x = i 

•• 

a px = 1 


Similarly for they 


and z directions: 


Yp 

u py 

a py 

Z P 

u pz 

a pz 


= r cos ( 0 + wt) 

= - r ( 0 + w) sin ( 0 + lot) + r cos ( 0 + cot) 

= r cos ( 9 + iot) - 2 r ( 9 + to) sin ( 9 + lot) 

• - r [ ( 9 + to)2 cos ( 9 + ^t) + 9 sin ( 9 + cot) J 
= r sin ( 9 + wt) 

= r ( 9 + to) cos ( 9 + tot) + r sin( 9 + tot) 

= r sin (.9 + cot) + 2 r ( 9 + 10 ) cos ( 9 + lot) 

+ r [ 0 cos ( 0 + lot) - ( 9 + to) 2 sin ( 9 + cot) ] 



(5) 


( 6 ) 
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The absolute components of the gas velocity may be expressed in 
terms of the gas relative velocity components. They become: 

u gx “ a gl 

Ugy ss - r to sin ( 3 + wt) - Ug 9 sin( 0 + wt) + u gr cos { 9 + wt) 
u gz - r a) cos ( 9 + wt) + Ug g cos ( 9 + wt) + Ug r sin( 9 + wt) (7) 

Substituting equations (4) - (7) into equation (3) and rearranging 
to solve for the particle relative acceleration components yields: 

1 = B (u gx - 1) (8) 

r = 3 (u gr - r) + r ( 9 + w) 2 (9) 

3 = £ (^ g Q -r 3 ) - 2 £ ( 9 + w) , 

r r 

Equations (8) - (10) constitute a coupled system of second order 
ordinary differential equations which govern the relative motion 
of a particle through a rotating blade row of a turbomachine. As 
pointed out in reference 7, the second term in the right hand 
side of equation (9) accounts for the centrifugal acceleration of 
the particle, while the second term in the right hand side of 
equation (10) represents the Coriolis component of acceleration. 
The path followed by a particle and its velocity along that path 
may be determined by numerically integrating equations (8) - (10) 
through the gas flow field. 

In the absence of reliable drag data the drag coefficient is 
frequently assumed to be given by the classical Stokes drag, law 
C D = 24/Re, where Re = D P g |ug-u p | /y , Although it is most con- 
venient to use Stokes law for the drag coefficient, it becomes 
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increasingly unrealistic for Reynolds numbers above approximately 
1.0 and in flows where rarefaction, compressibility and inertial 
effects are significant. The drag coefficient used in this study 
is given by Carlson and Hoglund (ref. 12) and corrects the Stokes 
law for these effects. Therefore, it is applicable throughout 
the continuum, slip and transition flow regimes. The expression 
is 

C D = 2A_ ip(Re, M) , (H) 

Re 

where 

i(j (Re , M) = (1 + .15 Re~ 687 ) (1 + exp (-. 427/M 4 • 65 -3/Re~ 88 )) 

1 + (M/Re) (3.82+ 1.28 exp (-1.25 Re/M) ) 

accounts for deviation from Stokes flow. The Mach number is 
based upon the absolute relative velocity of the gas with respect 
to the particles. 

PARAMETRIC STUDY 

The dynamic behavior of several gas-particle mixtures was 
determined by numerically solving equations (8) , (9)., and (10) 
with the viscous drag coefficient given by equation (11) . Several 
cases were examined by varying the diameter and mass density of 
the particles over ranges selected to correspond to values encoun- 
tered in most present-day LDV gas dynamic applications. The 
particle diameter varied from 0.5-4 microns and the mass density 
ranged from 1-4 gm/cc. The rotor gas flow remained constant 

throughout the study and was characterized by a mass flow rate 

3 

of 29.5 kg/sec (65 lbm/sec) , inlet mass density of 0.845 kg/m 
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(0.0527 lbm/ft 3 ) , inlet total pressure and temperature of 3.37 N/m 2 
(10.13 lbf/in 2 ) and 288 K(518.4R), respectively, and a specific 
heat ratio of 1.4. The rotor speed was constant at 16100 RPM. 
Particles were introduced into the flow system at initial radial 
positions extending from 60-94.3 percent of the radial span 
(. 19 8^r 4 *. 241 m) . The relative gas Mach numbers at the point of 
particle injection ranged from 1.15 at the 60 percent span point 
to 1.31 at the 94.3 percent span point. The particles were 
assumed to possess an initial velocity equal to 99 percent of the 
corresponding gas velocity and to have no angular deviation 
relative to the gas. Their velocities and trajectories were 
numerically determined as they progressed through the blade chan- 
nel. Analysis of a particlar particle was continued until either 
it passed through the rotor or until it entered the near-blade 
region. This is defined as the region between the finite differ- 
ence grid network and a blade surface. Since particles which 
entered this region were generally large and dense, and therefore 
unsuitable for LDV applications, their analysis was terminated 
immediately after entering the region. In an actual blade channel, 
particles which enter the near-blade region normally collide with 
the adjacent blade surface and rebound back into the flow stream 
with an increased velocity lag and angular deviation relative to 
the. gas. In practice, LDV electronics systems are generally de- 
signed to descriminate against such particles either because of 
their large size or because their velocities are significantly 
different from the mean velocity of the surrounding particles. 
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GAS FLOW RESULTS 


The three-dimensional gas flow field was determined with 
the numerical technique of reference 11 on blade-to-blade stream 
surfaces corresponding to 60,0, 69,2, 77.8, 86.2, and 94.3 percent 
radial span based upon the radius at the leading edge. Stream 
surfaces are approximated by conical surfaces of revolution with 
their apexes oh the compressor axis of rotation. They are char- 
acterized by radii at the leading and trailing edges of the blade 
row respectively, and by their cone angle. Table I summarizes the 
stream surfaces examined in this study. 

The flow fields on the 94.3, 77.8, and 60.0 percent span 
stream surfaces are illustrated in figures 3, 4, and 5 by relative 
intrablade Mach number vectors. They are located at each finite 
difference grid point within the solution region. The regions 
are presented in terms of nondimensional 


TABLE I 

GAS FLOW SOLUTION REGIONS 


mm irsm 

Radius at 
Leading Edge, m 

Radius at 
Trailing Edge, m 

Cone 

Angle, Deg. 

Percent 

Span 

RP2 

.242 

.238 

-8.227 

94.3 

RP3 

.232 

.229 

-6.120 

86.2 

RP4 

.222 

.220 

-3.946 

77.8 

RP5 

.211 

.210 

-1.705 

69.2 

RP6 

.200 

.202 

.642 

60.0 


meridional and tangential coordinates respectively. The axial 
and radial velocity components at each- grid point are determined 
from the mer.idonal component and the stream surface cone angle. 
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The gas velocity at any axial, radial and tangential position 
within the three-dimensional rotor blade channel is determined 
with an interpolation routine based upon flow field information 
on the stream surfaces of Table I. Table II summarizes average 
flow field information across the entrance and exit of the sol- 
ution regions. 


TABLE II 



Figures 3-5 provide justification for the following con- 
clusions regarding the three-dimensional gas flow field relative 
to the blade: 

a) the intrablade velocity increases with radial span, 

b) the inlet flow field is supersonic and the outlet flow 
field is generally subsonic, 

c) there is generally a smooth turning of the flow through 
the blade channel, 

d) there is a deviation in the smooth transition of flow angle 
near the leading edge of the pressure surface and the trail- 
ing edge of the suction surface. 
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e) On each stream surface the velocity decreases as the 
tangential position increases at constant meridional 
position from the suction surface to the pressure surface, 

f) the relative Mach number adjacent to the leading edge of 
the suction surface varies with radius from 1.28 to 1.69, 

. and the relative Mach number adjacent to the trailing edge 
of the pressure surface varies from 0.3465 to 0.687,. 

PARTICLE FLOW RESULTS 

The results of numerically integrating eqns. (8) - (10) 
through the rotating blade channel are summarized in figures 
6-20. Figures 6-8 present several particle trajectories 
projected onto the axial- tangential coordinate surface. The 
trajectories are shown relative to the blade, and the average 
radial coordinate of the trajectories is 0.220 m (.721 ft.). It 
will be subsequently shown that these trajectories lie essen- 
tially in the 77.8 percent radial span (RP4) gas stream surface. 
Figures 9 and 10 show several trajectories plotted on the axial- 
radial coordinate plane. 

Figures 6 and 9 present a family of trajectories of 1 micron 
diameter particles with a mass density of 1 gm/cc. The trajec- 
tories shown in figure 6 are typical of those obtained throughout 
the study and represent the approximate range of trajectories 
examined which did not intersect either the upper or lower near- 
blade region. The trajectories in figure 9 closely parallel the 
intersection of the gas stream surfaces of Table I with the axial- 
radial plane. (shown as dashed lines) . Therefore the particles move 
essentially on gas stream surfaces. 
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The influence of particle diameter and mass density is 
shown in figures 7, 8, and 10. Figures 7 and 8 indicate that an 
increase in either parameter causes an extension or lagging of 
the trajectory on the axial-tangential surface and this results 
in a corresponding increase in the velocity lag and angular devi- 
ation between the particles and the gas. While the results shown 
in Figures 7 and 8 are for a particle mass density and diameter 
of 1 gm/cc and 1 micron respectively, the results are typical of 
all other cases examined. It's seen that the curved particle 
path on the 1 - 9- surface is more sensitive to an increase in 
particle size than to a corresponding increase in mass density. 
This is justified by noting from equation (2) that if Stokes 
flow prevailed (<MRe, M)=l), the particle diameter and mass den- 
sity would appear to the second and first power respectively. 
Under these conditions a doubling of the diameter would have the 
same effect as a 4 to 1 increase in mass density. While Stokes 
flow does not always occur in this study, the above anology re- 
mains essentially valid since the drag expression given by equa- 
tion (11) is based upon the Stokes drag law. 

Figure 10 illustrates the superposition on the axial-radial 
plane of several particle trajectories resulting from particles 
of differing diameters and mass densities. As in the case of 
figure 9, these particle paths lie essentially on gas stream 
surfaces. The results clearly indicate that within the range of 
particle diameters and densities examined, the radial variation 
of a particle passing through a rotating blade channel is 
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independent of the size and. density of the particle. An examina- 
tion of the magnitude of the centrifugal acceleration term in 
equation (9) shows it to be small relative to the other term, 
thereby confirming .this result. 

The tangential variation of particle-to-gas velocity ratio 
and particle angular deviation at three axial positions is illus- 
trated in figures 11-14 as a function of particle diameter and 
average radial position. The abscissa is the normalized tan- 
gential position and is 0 at the suction surface and 1 at the 
pressure surface. The velocity ratio is defined as the magnitude 
ratio of the particle velocity to gas velocity V/U where V and U 
are the absolute particle and gas velocities respectively. Radial 
velocity components are not included in the determination of V 
and U since an LDV system normally views the blade channel from 
the radial direction, and is therefore not able to measure the 
radial component. Particle angular deviation * is defined as 
(j> = | Op - ctgl where the flow angles a p and Og are based on the 
absolute axial and tangential velocity components. These param- 
eters characterize the dynamic behavior of a particle at a partic- 
ular point in the fluid flow. Particle velocity lag is defined 
as | u — v|/U and is a related alternate parameter. Such a param- 
eter does not however indicate whether the particles are moving 
faster or slower than the surrounding gas. LDV applications typ- 
ically require that .99 « V/U * 1.01 and <j> '< .5 degree. The tan- 
gential profiles of V/U and 4> do not extend the full tangential 
. width of the stream channel because of the intersection of the 
outer trajectories with the near-blade region and the channeling. 
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(contraction) effect of trajectories with increasing axial stream- 
length and particle size. 

Figures 11-14 provide justification for the following con- 
clusions regarding particle velocity lag and angular deviation: 

a) particle velocity lag decreases (V/U -*■ 1) when the particle 
diameter decreases, when the radial position increases, or 
when the middle of the blade channel is tangentially ap- 
proached from either blade surface. 

b) particle velocity lag and angular deviation are generally 
worst near the pressure surface at the entrance to the 
blade channel. Angular deviation is also large near the 
suction surface for axial positions greater than approx- 
imately 1/L = 0, 

c) particle angular deviation becomes increasingly smaller 
as the particle size decreases and as the middle of the 
blade channel is tangentially approached from either blade 
surface, 

d) minimum angular deviation and velocity lag occur generally 
in the central portion of the blade channel where 1/L = 0 
and 0.4 ^9^ 0,6. 

An additional result observed, although not illustrated 
in figures 11-14, was that doubling the particle size causes ap- 
proximately the same change in V/U and $ as a 4 to 1 increase in 
mass density. Therefore, the particle path would remain essentially 
unchanged if the particle diameter was halved while at the same 
time increasing its mass density by a factor of 4. This result 
was also noted and graphically demonstrated in reference 9. 
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Figures 15-20 illustrate the tracking capability of particles 
in blade channels at three different average radial spans. The 
particle diameters varied from 0.5-2 microns and the mass density 
remained constant at 1 gm/cc. Tracking is graphically shown by 
lines of constant velocity ratio V/U and angular deviation <j> . 

The steepest gradients in V/U occur near the leading edge of the 
blade channel and in the region immediately before the leading 
edge. The high gradients in these regions are attributed to the 
sudden stream channel contraction and a "smeared" bow shock in 
front of the leading edge of the suction blade surface. 

The areas of greatest velocity lag occur generally in several 
well defined regions or pockets. One region is adjacent to the 
leading edge and another is just prior to the trailing edge of 
the suction surface. These are high blade-curvature regions and 
the particles are moving slower than the gas (V/U < 1) and are 
attempting to catch up. Two other high-lag areas are located 
immediately in front of the pressure surface leading edge and in 
a narrow oblique area in front of the blade channel. Particles 
in these respective areas are actually moving faster than the gas 
(V/U > 1) because of the sudden gas deceleration near the leading 
edge of the pressure surface, and deceleration across the "smeared" 
bow shock. 

Minimum velocity lag generally occurs near the entrance to 
the flow region where the gas-particle mixture was introduced and 
in the region beyond the blade channel. A broad area of good 
particle tracking is also seen to lie in the upper central portion 
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of the blade channel near the pressure surface. Figures 15-17 
indicate that velocity particle lags of 0.5 and 1 micron diameter 
particles are generally one percent or less in most regions of 
the blade channel. Two micron diameter particles are generally 
characterized by velocity lags of two percent or greater. These 
results remain essentially valid regardless of the radial position. 

Figures 18-20 indicate that the regions of greatest angular 
deviation occur near the leading edge of both blade surfaces and 
near the trailing edge of the suction surface. Minimum angular 
deviation occurs near the entrance to the flow region and in the 
central portion of the blade channel. A generalized conclusion 
regarding particle-gas angular deviation for most locations within 
the blade channel is that the angular deviation is 0.5 and 1 degree 
or less, respectively, for 0.5 and 1 micron diameter particles, 
and greater than 1 degree for 2 micron particles. 

CONCLUSIONS 

The results of this investigation have provided information 
regarding the dynamic behavior of micron size particles moving 
in the three-dimensional flow field of a rotating transonic axial- 
flow air compressor rotor. The results indicate that particle 
velocity lag and angular deviation are increased whenever particle 
size or mass density is increased and that these parameters are 
more sensitive to a change in particle size than to a corresponding 
change in mass density. A 2 to 1 increase in particle size was 
found to have essentially the same effect on particle tracking as 
a 4 to 1 increase in particle mass density. 
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Plotting of particle trajectories on the axial-radial plane 
show that the particles move essentially on gas stream surfaces 
and that their radial position within the blade channel is inde- 
pendent of particle size and mass density. Velocity lag and 
angular deviation are greatest in several regions near the leading 
and trailing edges of the blade channel. Favorable lag and de- 
viation occur near the entrance and exit of the solution region 
and in the central portion of the blade channel near the pressure 
surface. While particle and gas velocities are functions of 
radial position, particle tracking is not sensitive to the radial 
coordinate. 

Results indicate that LDV applications employing 1 gm/cc 
tracer particles with diameters greater than 1 micron would expe- 
rience velocity and angular deviations generally greater than 1 
percent and 1 degree respectively. While these conclusions are 
based upon a specific transonic flow through a rotating compressor 
rotor with a given blade channel geometry and a fixed speed, it 
is felt they are also approximately valid for other high speed 
flows through similar compressor rotor configurations. These 
generalized conclusions should not be applied to regions immedi- 
ately behind strong shock waves. Since shocks were assumed to 
be "smeared" in this study, there were no discontinuities in the 
gas velocity across the shocks and therefore, no high-lag relax- 
ation zones behind them. Reference 5 showed that a relaxation 
zone of finite length exists behind a shock wave within which 
momentum exchange causes a gradual reestablishment of kinematic 
equilibrium between the phases. The presence of such a relaxation 


22 



zone would cause the conclusions in this study to be invalid in 
this region. 

It should be emphasized that this study did not investigate 
particle flow in the boundary layers of the blade channel nor 
did it examine particle collision phenomena with the blade surfaces. 
Therefore, it is conceivable that there may be particles within 
the blade channel which have previously collided with a blade sur- 
face and possess a velocity and flow angle significantly differ- 
ent from the surrounding undeflected particles. However, refer- 
ence 9 showed that unless the particles are very large and dense, 
the dynamic behavior of blade impacted particles is essentially 
indistinguishable from the behavior of the surrounding undeviated 
particles. 
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APPENDIX-SYMBOLS 


a 

C D 

D 

L 

m 

M 

P 

Re 

t 

u 

U 


V 


x, y, z 
i, r, e 


. . . 2 
acceleration, m/sec 

drag coefficient 

particle diameter, microns 

blade chord, meters 

meridional coordinate, meters 

Mach number 

pressure, N/m* 

Reynolds number, p g D|u g " u pl/^' 
time , sec 
velocity, m/sec 

absolute gas velocity, composed of axial and tangential 
components only, m/sec 

absolute particle velocity composed of axial and tan- 
gential components only, m/sec 
fixed coordinate axes 

axial, radial, and tangential coordinate axes of blade- 
fixed coordinate system or components of the relative 


particle position vector 

1, r, 9 particle velocity components relative to the blade 

1, *r , 9 particle acceleration components relative to the blade 


a angle formed by absolute axial and tangential velocity 

components , deg . 

o -1 

3 flow parameter, 18 y t ( Re * M) /D pp, Sec 

y ratio of specific heats 

9 tangential coordinate , rad 

u gas viscosity, N-sec/m 
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p 


X 


^(Re, M) 

Gd 


density, kg/m 3 

angle between meridional streamline and axis of rota 
tion, deg. 

absolute angular deviation between the particles and 
gas, | Op - 0g| , deg. 

Reynolds -Mach number dependent function 
rotor angular velocity, rev/min 


Subscripts 


g 

P 

x, y, z 

1. r, 0 


gas 

particle 

absolute components 

axial, radial and tangential relative components 
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FIGURE 1 


Compressor Rotor 
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FIGURE’ 2 Blade-To-Blade Surface of Revolution 
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FIGURE 3 Gas Flow Field at 94,3 Percent Span 
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FIGURE 4 Gas Flow Field at 77,8 Percent Span 
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FIGURE 6 Typical Particle Trajectories 
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FIGURE 7 Particle Trajectories as a Function 
of Particle Diameter 
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FIGURE 8 Particle Trajectories as a Function 
of Particle Rass Density 
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FIGURE 9 Particle Trajectories on the Axial-Radial Plane 
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FIGURE 10 Particle Trajectories on the Axial-Radial Plane 
as a Function of Particle Properties 






FIGURE 11 Tangential Variation of Particle-To-Gas 
Velocity Ratio as a Function of Particle 
Diameter, r = - 220 m 

























FIGURE 14 Tangential Variation of Angular 
Deviation as a Function of Radial 
Position, d = x Micron 









FIGURE 15a Contours of Constant Particle-To-Gas Velocity 

RATIO, Pp = 1 gm/cc, r = .238 m 
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FIGURE 15b Contours of Constant Particle-To-Gas Velocity 

RATIO, Pp = 1 gm/cc, r = .238 m 
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FIGURE 15c Contours of Constant Particle-To-Gas Velocity 

RATIO, Pp = 1 gm/cc, r = .238 m 
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FIGURE 16a Contours of Constant Particle-To-Gas Velocity 
Ratio , p p * 1 9 m / cc ' r = .220 m 
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FIGURE 16b Contours of Constant Particle-To-Gas Velocity 

RATIO, pp = 1 gm/cc, r = .220 m 
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FIGURE 16c Contours of Constant Particle-To-Gas Velocity 

RATIO, P p ■ 1 gm/cc, r «* . 220 m 
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FIGURE 17a Contours of Constant Particle-To-Gas Velocity 

RATIO, Pp = 1 gm/cc, r = .200 m 
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FIGURE 17b Contours of Constant Particle-To-Gas Velocity 

RATIO, P p = 1 gm/cc, r = .200 m 
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FIGURE 17c Contours of Constant Parti cle-To-Gas Velocity 

RATIO, Pp = 1 gm/cc, r = .200 m 


44c 



THETB/THETFUREF) 



(BXIRL STRERMLENGTH.J/L 


FIGURE 18a Contours of Constant Angular Deviation, 

Pp = 1 gm/cc , r = ,238 m 
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FIGURE 18b Contours of Constant Angular Deviation, 

Pp = 1 gm/cc, r --.238 m 
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FIGURE 18c Contours of Constant Angular Deviation, 

Pp = 1 gm/cc f r = .238 m 
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FIGURE 19a Contours of Constant Angular Deviation., 

Pp ss 1 gm/cc, r = .220 iin 
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FIGURE 19b Contours of Constant Angular Deviation, 

Pp = 1 gm/cc, r - .220 m 
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FIGURE 19c Contours of Constant Angular Deviation, 

Pp '= 1 gm/cc, r = .220 m 
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FIGURE 20a Contours of Constant Angular Deviation, 

Pp - 1 gm/cc, r = .200 m 


47a 



FIGURE 20b Contours of Constant Angular Deviation, 

Pp = .1 gm/cc, r = ,200 m 
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FIGURE 20c Contours of Constant Angular Deviation, 

pp = 1 gm/cc, r = .200 m 
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